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Soft tetragonal distortions in ferromagnetic Ni2MnGa and related materials from first 
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A detailed examination of the energy landscape, density of states and magnetic moment of tetrago- 
nally distorted ferromagnetic Ni2MnGa was performed using first-principles local-spin-density (LSD) 
pseudopotential calculations, varying V as well as c/a. The energy of tetragonal Ni2MnGa is found 
to be nearly constant for values of c/a over a wide range, with shallow minima near c/a = 1 and 
1.08 in addition to that near 1.2. This flat energy surface is consistent with the wide range of 
observed values of c/a. It also explains the observation of pseudomorphic growth of Ni2MnGa on 
GaAs, despite a nominal 3% lattice mismatch. The related materials Ni2MnAl, Ni2MnIn and fer- 
romagnetic NiMn were examined for similar behavior, but all are seen to have a single well-defined 
minimum at c/a near 1, consistent with available experimental information. For NiMn, the ground 
state antiferromagnetic ordering and structural parameters are correctly predicted within the LSD 
approach. 
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In the design of mechanical actuator devices, both on 
macroscopic scales and in thin- film based MEMS, there 
is continuing interest in identifying and optimizing new 
high-performance materials. Materials which exhibit a 
martensitic transformation and associated shape memory 
effect have been shown to be quite useful, though their 
high-frequency applications are limited by the slow rate 
of the martensitic transformation and poor energy con- 
version. It has been proposed^ that this deficiency could 
be addressed by using shape memory materials which are 
ferromagnetic, and using appJifid magnetic fields to con- 
trol the mechanical responscnQ. 

Therefore, it is of particular interest to identify can- 
didate ferromagnetic martensites and to understand the 
transition mechanism, especially the coupling between 
magnetic and structural degrees of freedom. The rather 
small number of known systems includes Fe-Pt, Co- 
Ni, Fe-Mn-C, Fe-Ni-C, Fe-Ni-Co-Ti and Fe-Ni, with the 
Heusler structure compound NipMnGa being the most 
thoroughly investigated to dateail2l. 

In this paper, we investigate the structural energetics 
of Ni2MnGa and related materials from first principles, 
with a particular focus on the tetragonal distortion asso- 
ciated with the martensitic transition. For comparison, 
we examine Ni2MnAl and Ni2MnIn (not observed to be 
martensitic), and NiMn, related by replacing Ga with 
an additional Mn (exhibits a martepsitic transition to an 
antiferromagnetic tetragonal phasdHl). While our calcu- 
lations are for bulk crystals, they also allow us to predict 
the effects of epitaxial stress, which often is the dominant 
factor in determining the structure and properties of thin 
films. 

We use_tjhe self-consistent pseudopotential plane wave 
approachEj within the local spin density approximation 
(LSD A). We use Troullier-Martinai3 pscudopotentials 
and exchange and correlation potential in the Perdew- 



TABLE I. Reference configurations and cut-off radii (a.u.' 
used to construct the pscudopotentials. 
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Ni 3d''4s^4p'-' 


2.2 


2.2 


2.2 


Mn 3d"4s''4p" 


1.9 


2.6 


2.0 


Al 4sMp"4d" 


2.3 


2.3 


2.3 


Ga 4s^4pMd" 


3.0 


2.6 


3.0 


In AsMp'M" 


3.0 


3.0 


3.0 



WangEJ form. The nonlinear core correction schemeta 
was used in the construction of the pscudopotentials. 
The pseudopotential cut-off radii are summarized i^i-i Ta- 
ble I. The local components of the pseudopotentialsEEl are 
chosen to be 1=0 for Ni, Mn, Al and In and 1=1 for Ga. 
The energy plane-wave cut-off is 70 Ry. The unit cells 
consist of one formula unit for Ni2MnX (X = Al, Ga, 
and In) and two formula units for NiMn. The Brillouin 
zone ia.sampled by a 10x10x10 k-point Monkhorst-Pack 
(MP)EZl mesh with zero shift. To be able to calculate the 
magnetic moment accurately, we do not apply tempera- 
ture smearing for the k-point integration. In NiMn, we 
studied two antiferromagnetic phases: AF-I, with the Mn 
atoms with parallel moments located on alternating (1 1 
0) planes; and AF-II, with parallel moments on alternat- 
ing (1 1 1) planes (in the cubic reference system)lla. The 
self-consistent calculations tend to converge to a local fer- 
romagnetic minimum. To make the system converge to 
the antiferromagnetic state, approximately 0.2 eV/atom 
lower in energy, we constrain the total magnetic moment 
to zero. 

For the high temperature cubic structures, our calcu- 
lations yield the lattice constants, bulk moduli and mag- 



TABLE II. Lattice constants, bulk moduli and mag- 
netic moments of Ni2MnX compounds (L2i structure) and 
FM-NiMn (B2 structure) compared to experiment (in paren- 
theses) and previous theoretical calculations (square brack- 
ets), where available. 
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FM-NiMn 




5.56 (5.63") 


NiaMnAl 


10.93 


(11.01") [10.98"] 


NiaMnGa 


10.91 


(11.01") [10.95"] 


NiaMnIn 


11.32 


(11.47") [11.43=] 


B (GPa) 


FM-NiMn 




155 


NiaMnAl 




167 [163"] 


NiaMnGa 




170 (146") [156"] 


NiaMnIn 




138 


H{^B/Mn) 


FM-NiMn 




4.42 [3.83] 


NiaMnAl 


4.22 (4.19") [4.03"] 


NiaMnGa 


4.22 (4.17") [4.09"] 


NiaMnIn 


4.31 (4.4C 


)") [3.91"] [4.16-'] 
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netic moments given in Table IT For NiMn, we compare 
calculated results for a hypothetical ferromagnetic cu- 
bic phase with measurements in the paramagnetic cubic 
phase. The results compare well with experiment and 
with previous first-principles studies. 

Next, we present the results of total energy calcula- 
tions for uniform tetragonal distortions of these cubic 
structures, and discuss them in the context of observed 
low-symmetry phases. At least three distinct ferromag- 
netic low-symmetry phases have been experimentally ob- 
served in bulk Ni2MnGakl, two tetragonal (/3' and /?'") 
and one orthorhombic (/3"). The (3' phase, obtained by 
cooling of the high-temperature Heusler (L2i) phase be- 
low Tm ~ 200 K, has a tetragonaLstructure with c/a = 
0.94 (a=5.920 A and c=5.566 A).eI The relative volume 
change across the martensitic transition is only 1%. In 
addition, there is a shufHe of the (110) planes along 
[110]. This has incommensurate periodicity represented 
by a wavevector of 7r/a[0. 43, 0.43,0], corresponding to ap- 
proximately five interplanar distancea3'B. The transition 
is preceded by a cubic premartensitic {L2pm) phase be- 
low Tpj\/ w 260 K, characterized by softening and con- 
densation of transverse acoustic phonons 7r/a[0. 33, 0.33,0] 
with displacements similar in type to the shuffle distor- 
tion below Tm- Under uniaxial compression along [1 1 0] 
the f3' phase transforms into the orthorhombic f3" phase 
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FIG. 1. Total energy as a function of c/a in NiaMnGa. 
Each curve corresponds to a constant volume V per unit cell, 
measured relative to V/feusier ~ 1300 a.u."^ The curves are 
separated by intervals of AV/V = 1.5%. Each curve repre- 
sents an interpolation through 26 points, shown explicitly by 
filled circles for V = V HE-a.siar- 



(a = 5.54 A, b = 5.78 A and c = 6.12 A)Q. X-ray analy- 
sis shows that /?" martensite is also modulated with [lIO] 
displacements of the (1 1 0) planes with a periodicity of 
seven (110) interplanar distances. The third martensite 
phase /3"', obtained by further compression along [110], 
is metastable with respect to the removal of the applied 
stressQ. The /?'" structure is a uniform tetragonal dis- 
tortion of the cubic Heusler structure (a = 6.44 A and c 
= 5.52 A) with c/a = 1.18. The relationships between 
these phases can be summarized as follows: 



L2i 
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The calculated energy of Ni2MnGa as a function of c/a 
for V — V Heusler IS showu ill Figurc |l|. The region of c/a 
from 0.9 to 1.3 is examined in considerable detail, with 
26 separate calculations. Based on the experimental ob- 
servations of bulk crystal structures, one would expect to 
find a local minimum in the energy as a function of c/a 
near 1.18 and perhaps another near 0.94. There is a shal- 
low low-energy local minimum near c/a =1.19 with total 
energy 2 meV/atom higher than the L2i structure. In 
the previous calculations of Ayuela et.al., this minimum 
is at c/a=1.16 and is slightly lower in energy than the 
L2i structure. The minimum can readily be associated 
with the pure tetragonal f3"' phase observed by Kokorin 
et.al. with c/a=1.18Ll. In contrast, there is no local min- 
imum or any discernible anomaly at c/a — 0.94, with a 
smooth decrease to the minimum at the L2i structure. 
However, our fine-scale exploration of the energy surface 
does reveal that the entire energy surface for 0.95 < c/a 
< 1.25 is remarkably flat, with the total energy varying 
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FIG. 2. Total energy as a function of c/a ratio of several 
related compounds, compared with Ni2MnGa. The curve for 
each compound is computed at a constant volume equal to the 
theoretical V//eus!er given in Table II. Each curve represents 
an interpolation through at least 20 points in the range shown. 
The energies are given relative to the energy at c/a = 1, with 
the curves slightly offset for clarity. 



less than 2.5 meV/atom. This has significant physical 
consequences for thin films which will be discussed fur- 
ther below. 

First, we consider the reasons for the absence of a local 
minimum corresponding to the ground state /3' phase. It 
has been suggested that the tetragonal distortion could 
be stabilized by a change in volumeo. To investigate this 
hypothesis, we computed the energy as a function of c/a 
for varying volumes, with results shown in Figure |lj. We 
sampled the volumes in 1.5% increments of /S.V/V . Even 
with the possibility of volume change, the L2i structure 
remains the total energy minimum, with no local mini- 
mum corresponding to the (3' structure. We remark that 
the minimum near c/a = 1.19 is lowered relative to the 
cubic structure as V increases, although it never actu- 
ally becomes more favorable. We have also calculated 
the total energy of a pure orthorhombic structure at the 
experimental lattice constants, which give c/a=l.ll and 
b/a=1.04.IZl The calculated energy of this orthorhombic 
phase is 4 meV per atom higher than the L2i phase, con- 
sistent with previous calculations. El We conclude that the 
observed shuffle distortions are crucial to stabilizing the 
observed /?' and /3" phases. 

A number of related materials were examined for com- 
parison. Martensitic transitions have been observed 
in the NijMn-Al system, but not in stoichiometric 
Ni2MnAlL2ltJ. Similarly, there are no reports of low- 
symmetry phases of Ni2MnIn. NiMn exhibits a marten- 
sitic transition from a cubic paramagnetic phase at Tm 
= 943 K to a tetragonal phase, wiih c/a = 1.33 which, 
however, is antzferromagnetictijEj. 

For Ni2MnAl, Ni2MnIn, and ferromagnetic NiMn, the 



FIG. 3. Total energy of NiMn (FM, AF-I, AF-II phases) as 
a function of c/a at constant volume. For FM, AF-I and AF-II 
phases, w^ji-eJiose the volume of the observed experimental 
structureBN, a=3.74 A and c=3.52 A. 



energy as a function of c/a for V — Yneusier is shown in 
Figure 0. All three compounds are seen to have a sin- 
gle well-defined minimum at c/a near 1, consistent with 
available experimental information. In the previous cal- 
culations of Ayuela et.al., there was another metastable 
minimum for Ni2MnAl at c/a=1.22. We find that for 
Ni2MnAl the shape of this curve is very sensitive to the 
k-point sampling density, with a local minimum near 
c/a=1.2 evolving into a slope "discontinuity" for grids 
denser than 8x8x8. 

When constrained to be ferromagnetic, NiMn has a 
tetragonal distortion energy curve similar in shape to 
that of Ni2MnIn and Ni2MnAl, though with a larger 
stiffness. The substitution of a magnetic atom in the 
X site has, however, significant consequences for the 
magnetic order, stabilizing a type-I antiferromagnetically 
ordered structure (AF-I) by 190 meV over the lowest- 
energy ferromagnetic structure (FM). This is in good 



agreement with the_.value of 200 meV obtained in a 
previous calculatiorJla, performed with lattice constants 
fixed at experimental values. As shown in Figure 0, to- 
tal energy minimization for the AF-I NiMn structure 
leads to a pimdicted value of c/a — 1.33, consistent with 
experiment 0'E£l At this value a pseudogap appears at the 
Fermi level, as shown in Figure 11. As the tetragonal dis- 
tortion deviates from c/a = 1.33 the pseudogap becomes 
less pronounced and finally Ni d-states fill the gap. The 
formation of a pseudogap and the associated enhanced 
stability does not appear for a distinct antiferromagnetic 
ordering (AF-II), which is higher in energy by 70 meV 
and has a much shallower minimum at c/a = 1.24. 

The densities of states for cubic Ni2MnX (X — Ga, Al 
and In) are shown in Figure |5[ For all three compounds, 
these are quite similar, especially near the Fermi level. 
The total magnetic moment of Ni2MnGa L2i structure is 





FIG. 4. Density of states for (a) the AF-I and (b) the AF-II 
phases of NiMn for varying c/a at the same volume as in 
Figure |3[ The sohd hne corresponds to the minimum energy 
states for both AF-I and AF-II phases. 




FIG. 5. Density of states in Ni2MnX compounds (X — Al, 
Ga, In). The atomic character (obtained by projecting the 
wave functions onto individual atoms) is indicated by appro- 
priate labels. 




FIG. 6. Total magnetic moments of Ni2MnX (X = Al, Ga, 
In) and FM-NiMn as a function of c/a, corresponding to the 
total energy calculations in Figure 0. 



4.228 ^lB, which decomposes into 0.301/3. 700/-0.075 hb 
on Ni/Mn/Ga atoms. This agrees well with the previous 
calculations of References g and 10 . 



The magnetic moment as a function of c/a, shown in 
Figure |6|, is again similar for all three compounds. It 
shows a sharp minimum at c/a near 1 and and a local 
maximum at c/a ~ 0.94. From the atomic decomposition 
of this curve for Ni2MnGa in Figure M, we see that the 
features can be directly attributed to the Ni contribu- 
tion, most likely associated with a small Ni-derived peak 
in the minority density of states located just below the 
Fermi level. It is interesting to note that the position of 
this maximum coincides with the value of c/a for the /?' 
martensitic structure. 

The softness of the tetragonal distortion in Ni2MnGa 
has significant implications for the growth and properties 
of thin films. Efforts to grow single-crystal thin films on a 
GaAs substrate (wi|th|-3|.l% lattice mismatch) have been 
carried out by MBEc3E3. Using a Sco.3Ero.7As interlayer, 
a tetragonal-structure film with lattice parameters a = 
0.565 nm and c = 0.612 nm was obtained. This is ex- 
actly the structure we expect based on the total energy 
curves in Fig. n^. As the in-plane a is reduced to 0.565 nm 
to match the substrate, the c axis increases to maintain 
the volume of the cubic Heusler phase, which minimizes 
the total energy over the whole range of c/a considered. 
The energy of the resulting structure, at c/a = 1.08, is 
less than 2 meV/atom higher than the cubic structure, 
explaining the observation that the thickness of this pseu- 
domorphic film, 300 A, far exceeds the critical thickness 
expected for the given lattice mismatch. While the large 
epitaxial stress should affect the structure and properties 
significantly, if the film can be released it should exhibit 
a martensitic transition and shgipe-memory properties. 

In another recent experimentnZl a 450 A thick film of 
Ni2MnGa was grown by MBE on a GaAs substrate using 
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a NiGa interlayer. The structure of this film was reported 
to be tetragonal, with c/a = 1.05 (a=5.79 A and c=6.07 
A). In this case, the lattice constant of Ni2MnGa is a 
close (0.6%) match to the interlayer, leading to the ex- 
pectation that the structure of the film should be cubic. 
Instead, the c parameter is expanded, resulting in a vol- 
ume 3% larger than the volume of the L2i structure. One 
explanation might be a slight change in the stoichiome- 
try of the film. Modelling the effect of this shift as an 
negative applied hydrostatic pressure, we now focus on 
the ISy/V — -1-3% curve in ^ The flatness of this curve 
is similar to that for V = ^Heusier- Thus, although the 
expanded system is formally mismatched to the NiGa in- 
terlayer, the constraint a = 5.79 A is associated with an 
energy cost of less than 2 meV/atom, which would allow 
pseudomorphic growth of this tetragonal film far beyond 
the critical thickness. 

This study highlights the fact that the energy surface 
of Ni2MnGa is far from simple. Studies of additional 
instabilities of the cubic Heusler structure via density- 
functional perturbation theory are in progress. In par- 
ticular, the energies of shuffle distortions as a function of 
wavevector can be obtained with this method for a bet- 
ter understanding of the observed j3' ground state struc- 
ture. This information wifl also allow us to model the 
properties of the high-temperature cubic Heusler phase 
within an effective Hamiltonian approach. We expect 
that these properties reflect the presence of large local 
distortions around the average cubic structure and there 
should be interesting differences from the properties of 
other Heusler phases whose ground state structure is cu- 
bic. 

In conclusion, we performed first principles calcula- 
tions of the total energy as a function of tetragonal dis- 
tortion (c/a) and volume V for Ni2MnGa and related 
compounds Ni2MnAl, Ni2MnIn, and NiMn. The total 
energy of Ni2MnGa at constant volume is remarkably 
flat in the range 0.95<c/a<1.25, varying less than 2.5 
meV/atom. This provides an explanation of the sur- 
prisingly high compliance of single-crystal fllms grown 
by MBE. In contrast, the related materials Ni2MnAl, 
Ni2MnIn, and ferromagnetic NiMn have a single well- 
defined minimum near c/a = 1. The densities of states 
and c/a dependence of the magnetic moments are quite 
similar in all compounds studied, leading us to conclude 
that the unique behavior of Ni2MnGa is the result of 
fine tuning. It may be possible to achieve this in other 
systems through appropriate alloying. 
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Palmstr0m, J. Dong, C. Bungaro, X. Huang and A. 
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FIG. 7. Total magnetic moment and contributions associ- 
ated with individual elements atoms in Ni2MnGa as a func- 
tion of c/a for V = V Heusler- Note the difference in the ver- 
tical scale between the first plot and subsequent plots. 
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